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PREFACE 


This  RAND  Memorandum  is  devoted  -to  the  examination  and  development 
of  new  statistical  techniques  for  analyzing  accelerated  life  test  data. 
It  is  a  part  of  RAND's  continuing  interest  in  the  theoretical  and 
practical  aspects  of  reliability. 

This  Memorandum  is  addressed  to  statisticians,  data  analysts, 
reliability  engineers,  and  others  interested  in  the  analysis  of 
accelerated  life  test  data. 

Professor  Richard  E.  Barlow,  University  of  California,  Berkeley, 
is  a  consultant  to  The  RAND  Corporation. 


SUMMARY 


This  Memorandum  examines  and  develops  new  techniques  for  analysing 
life  test  data.  It  obtains  estimates  for  the  life  distribution  in  the 
use  environment  based  on  data  from  both  the  use  and  the  accelerated 
environments.  The  techniques  require  only  physically  piausibl-  is- 
sumptions,  not  the  usual  ones  involving  specification  of  a  far  ::  y  of 
parametric  probability  distributions.  Procedures  are  given  for  testing 
those  assumptions  that  are  made. 
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1 .  INTRODUCTION 


It  is  common  practice  in  certain  life  testing  experiments  to 
subject  test  items  to  overstress  conditions.  This  may  take  many 
different  forms.  For  example,  one  might  subject  electronic  capacitors 
to  a  high  voltage,  ball  bearings  to  a  high  load,  or  a  mechanical  assembly 
to  a  strong  vibration.  The  purpose  is  to  shorten  the  time  to  failure. 

The  problem  is  to  predict  the  time  to  failure  in  the  norm.-  or  use 
environment  on  the  basis  of  such  accelerated  life  test  data.  •’See 
Winter,  Denison,  Hietala  and  Greene  (1964)  for  engineering  details  and 
bibliography. 1 

Let  Y  (X)  be  a  random  variable  with  distribution  G  (F)  vhere  Y  (X) 
denotes  the  time  tc  failure  under  normal  (accelerated)  conditions. 

Suppose  that  F  and  G  are  related  by  a  time  transformation  <y(t) ,  where 

F(t)  -  G|-<y(t)l 

so  that,  assuming  G  *  exists, 

O'(t)  **  G"lF(t)  . 

Many  authors  [e.g.,  Bessler,  Chemoff,  and  Marshall  (1962)]  assume 

or(t)  ■  o>t  , 

i.e.,  o(t)  is  a  scale  transformation  of  the  time  axis.  They  also 
assume  that  F  and  G  belong  to  specified  parametric  families-  usually 
the  exponential  distribution.  In  certain  cases  there  seems  o  be  some 
justification  for  assuming  that  the  form  of  the  acceleration  function 
(i.e.,  the  time  transformation)  is  known.  Examples  are: 
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1.  heaters  for  vacuum  tubes  [R.  L.  Guild  (1952)]; 

2.  ball  bearings  [g.  Lieblein  and  M.  Zelen  (1956)]; 

3.  paper  capacitors  [G.  J.  Levenbach  (1957)];  and 

4.  certain  transistors  [G.  A.  Dodson  and  B.  T.  Howard  (1961)]. 

The  novelty  of  the  problem  of  course  vanishes  when  or(t)  is  assumed 

known  and  the  problem  becomes  a  classical  life  testing  problem.  In 
this  study,  we  do  not  assume  that  or(t)  is  known--only  that  <*(t)  and 
the  life  distributions  satisfy  certain  geometric  restrictions  which 
we  believe  are  intuitively  acceptable  ami  reasonable  for  many  appli¬ 
cations.  However,  we  do  assume  that  some  sample  data  from  items  tested 
in  the  use  environment  are  available,  albeit  possibly  scanty. 

In  this  paper,  we  develop  least  squares  estimators  for  life  dis¬ 
tributions  based  on  the  sampling  distribution  assuming  that  the  failure 
rate  is  increasing  on  the  average.  Given  life  test  data  from  the  ac¬ 
celerated  and  unaccelerated  modes,  we  simultaneously  estimate  both 
distributions  assuming  only  that  the  life  distributions  have  increasing 
failure  rate  on  the  average  and  that  accelerated  items  tend  to  fail 
sooner  than  unaccelerated  items.  Applications  to  fatigue  data  and 
accelerated  life  teBts  on  duplex  capacitors  are  considered.  Procedures 
for  testing  the  underlying  assumptions  are  described. 
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2.  A  STATISTICAL  MODEL  FOR  ACCEI£  RATED  LIFE  TESTING 

We  assume  that  observations  Y  -  (Y,  ,  Y_  ,  Y  )  are  available  on 
items  in  the  unaccelerated  mode  having  distribution  G.  These  may  be 
results  from  a  truncated  or  censored  life  test.  It  is  essential,  how¬ 
ever  ,  that  we  have  some  failure  information  on  life  times  In  the  un¬ 
accelerated  mode.  In  the  accelerated  mode  it  will  normally  be  possible 
to  obtain  many  more  observations,  say  X  ■  (X^ ,  X£ ,  •••,  Xr)  from  a 
distribution  F  also  unknown,  where  again  these  may  result  from  a 
censored  or  truncated  test.  There  are  certain  natural  restrictions 
which  we  would  want  to  impose  on  the  time  transformation,  a,  relating 
F  and  G.  Since  acceleration  reduces  the  time  to  failure,  we  would 
expect  cr(t)  2-  t.  Clearly,  a  will  also  be  nondecreasing.  On  the  basis 
of  these  properties  alone  we  could  consider  the  problem  of  obtaining 
an  estimate  for  G  which  maximizes  the  joint  likelihood  and  satisfies 
the  restrictions.  This  problem  has  been  solved  by  Brunk,  Franck, 

Hanson,  and  Hogg  (1966).  However,  these  estimates  impose  no  re¬ 
strictions  on  the  failure  distributions.  Utilization  of  additional 
a  priori  Information  concerning  the  life  distributions  should  result 
in  Improved  estimates. 

Birnbaum,  Esary,  and  Marshall  (1966)  have  characterized  the  smallest 
class  of  failure  distributions  containing  the  exponential  distributions 
which  is  closed  under  the  formation  of  coherent  structures.  This  class 
is  precisely  the  class  of  distributions  with  increasing  failure  rate 
average  (IFRA) ;  i.e.,  if  the  failure  rate  r(t)  exists,  then—  f*  r(u)du 
is  nondecreasing  in  t.  More  generally,  a  distribution  F  is  IFRA  if 


0>  >n> 


-  l°g[  1_  -  nonc)ecreas  ing  in  t  2  0.  Obviously,  this 

the  class  of  distributions  with  nondecreasing  failure  rate 
the  exponential  distribution. 


includes 
at.d  hence 


STATEMENT  OF  THE  PROBLEM 

Given  ordered  observations  Y.  SY.fi  ...  fi  Y  from  a  distribution 

12  m 

G  (the  unaccelerated  life  distribution)  and  ordered  observations 

X,  fi  X„  fi  ...  fi  X  from  a  distribution  F  (the  accelerated  life  dis- 
12  n 

tribution) ,  we  calculate  the  empirical  distributions  G  and  F  corre- 
'  *  r  m  n 

A 

spondlng  to  G  and  F,  respectively.  We  wish  to  obtain  estimates  G^  and 
such  that 

1.  G  and  F  are  IFRA: 

m  n 

2.  G  (x)  fi  ?  (x)  for  all  x  (i.e.,  stochastic  ordering); 

m  n 

A  A 

3.  G  and  F  are  closest  to  G  and  F  ,  respectively,  in  a  least 

m  n  m  ns  r  J 

squares  sense,  which  will  be  made  precise  in  Sec.  5. 

will  then  be  our  estimate  for  the  life  distribution  in  the  un- 
m 

accelerated  mode  using  data  from  bo th  the  accelerated  and  unaccelerated 
environments . 

Since  in  applications  one  does  not  always  have  complete  life  test 


data,  we  also  consider  estimation  procedures  for  incomplete  data. 
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3.  LEAST  SQUARES  ESTIMATES  FOR  I FRA  DISTRIBUTIONS 

Assume  that  each  of  n  units  is  observed  over  some  or  all  of  its 
life.  Thus  unit  i  is  observed  until  it  either  fails  at  a  random  age 
T  or  has  attained  age  L^(is  lost  to  observation),  whichever  occurs 
first  (i  =  1,  2,  n) .  The  ,  called  limits  of  observation,  are 

constants  or  values  of  other  random  variables,  which  are  assumed  to  be 
independent  of  T^.  In  accelerated  life  testing  censored  and/or 
truncated  samples  are  especially  common,  so  that  it  is  important  to 
consider  such  data.  We  now  give  a  procedure  for  obtaining  an  estimate 
of  the  unknown  failure  distribution,  F,  assuming  that  -log  F(t)/t  is 
nondecreasing  on  [0,  •)  where  F(t)  »  1  -  F(t). 

3.1  PROCEDURE  FOR  OBTAINING  THE  LEAST  SQUARES  ESTIMATES 

(1)  Suppose  k  (1  s.  k  £  n)  failures  are  actually  observed.  Let 
s  Xj  s  ...  £  X^  denote  the  ordered  ages  at  failure;  for  convenience, 
define  X  *  0.  Let  n,  be  the  number  of  items  under  observation  just 
before  X  ^ ,  6j  be  the  number  of  failures  at  X j ,  and  be  the  number  of 
"losses"  in  [X^_^,  Xj)«  Then  3  nj  “  “  ^j+i*  Ordinarily,  6^ 

will  be  one,  but  may  be  more  than  one  if  the  sample  contains  tied  ob- 


J 


j  +  1.  Define 


servations.  If  there  are  no  losses  or  ties,  n 
the  product-limit  estimate  as 

1,  0  S  t  <  Xj 

Fn(t)  -  J  fi  (n j  -  SjjAy  X£  s  t  <  X1+1,  i  -  1,  E-l 
^  ^  ,t 


(3.1) 


where  E  =  | 

'  n 


°  ,  t  ^ 

k  if  k  <  n 
if  k  =  n  * 


Undefined  if  k  <  n. 
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This  estimate  coincides  with  the  usual  empirical  sampling  distribution 
estimate  if  there  are  no  losses  and  if  k  =  n.  [See  Kaplan  and  Meier 
(195P) . ] 

A 

(2)  We  obtain  a  least  squares  estimate,  F^(t),  for  the  failure 
distribution,  F,  requiring  the  estimate  to  have  the  IFRA  property 
assumed  for  the  distribution.  Let 


vy  ■  fn(xr>  -  vy>  • 

It  is  easy  to  verify  from  (3.2)  that  \n(Xp  will  be  nondecreasing 

in  i.  In  Sec.  5  we  show  that  this  estimate  is  closest  to  X  (t)  in  a 

n 

least  squares  sense.  Any  nondecreasing  function  defined  at  order 
statistics  as  in  (3.2)  is  obviously  a  permissible  estimate  of  X.  How¬ 
ever,  for  definiteness,  we  define 

0  ,  t  <  xx 

\i(t)  *  VV*  xi  5  t  <  Xi+1,  i  -  1,  ...,  E-l  (3.3) 

+  *  ,  t  *  xE. 

Our  estimate  for  the  distribution  function  will  then  be 
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Fn(t) 


t  <  X, 


expL  — X n (X A ) t '  ,  X  £  t  <  X  ,  i  =  1,  ....  E-l  (3.4) 


s  V 


3.2  PROCEDURE  FOR  OBTAINING  STOCHASTICALLY  ORDERED  ESTIMATES 

Given  observed  failures  Xj  s  X9  s  ...  s  ^  (k  s  n),  as  before, 
from  a  sample  of  size  n  from  a  distribution  F,  and  observed  failures 
Y^  £  Y?  £  . . .  £  Y^  ( r  £  m)  from  a  sample  of  size  m  from  a  distribution 
G,  wo  wish  to  simultaneously  estimate  F  and  G  assuming: 


i)  F  and  G  are  I FRA ; 
ii)  F(t)  £  G(t)  for  all  t  >  0. 


Let  F  and  G  denote  the  product  limit  estimates  for  F  and  G  as  pre- 
n  m  r 

viously  described  (cf.  Eq.  (3.1)).  Let  X^X^)  =  -log  Fn^i^Xi  and 

y  (Y.)  =  -log  G  (Y.)/Y,  denote  estimates  of  X(X  )  =  X.  and  of 
ra  j  °  m  j  j  i  l 

y(Y  )  E  Y j •  respectively.  We  wish  to  determine  the  X^  and  y^  which, 
subject  to  the  constraints  (3.6),  minimize 

j  h,  -  X„(X  )]2  F n{Xi1  +  I  [V  -  V.(Y  )]2  GjY  ).  (3.5) 

i*l  j«l  J  J  J 


The  minimization  must  be  performed  subject  to  the  constraints 


OSXj  i  s  ...  £  Xk, 

0  *  Y1  *  V2  *  ***  *  V 


(3.6) 


and 

A  A 

G(t)  :>  F  (t). 
n  n 
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To  understand  the  nature  of  the  stochastic  ordering  restriction 

it  will  be  convenient  to  relabel  the  Y^'s.  l*et  s  •  ••  s  Y^  <  ^ 

denote  the  first  block  of  Yj's--those  less  than  X^.  The  set  may  be 

empty.  Similarly  let  X.  .  £  Y.,  <  ...  s  Y. ,  <  X  oenote  the  block  of 

1—1  11  lj£  * 

Y , 's  in  the  interval  (X^  ^ ,  X^). 

Our  problem  is  to  minimize  (3.5)  subject  to  (3.6)  and 

V,,  £  X,  (i  *  1,  2,  ...,  k).  This  is  a  straightforward  quadratic 

iji  i 

programming  problem.  A  RAND  code  (RSQPF4)  is  available  for  solving  this 
problem.  In  some  applications  we  can  obtain  the  solution  to  our  problem 
in  a  manner  similar  to  that  in  which  we  obtained  tho  IFRA  estimate  pre¬ 
viously.  To  be  specific,  consider  all  possible  "interlacings"  of  the 

y.'s  among  the  X's.  For  each  such  linear  ordering  we  can  obtain  an 
,i  i 

explicit  solution  to  the  problem  of  minimizing  (3.5),  subject  to  the 
linear  ordering,  using  Theorem  1  of  Sec.  5.  Now  choose  that  solution 
corresponding  to  a  linear  ordered  constraint  set  which  has  minimum  sum 
of  squares.  (Unfortunately,  for  even  moderate  k  and  r,  the  number  of 
interlacings  that  must  be  considered  will  be  quite  large.) 


3.3  APPLICATION  TO  FATIGUE  DATA 

To  illustrate  the  ideas  discussed  above  we  consider  some  fatigue 
data  extracted  from  a  paper  by  R.  P,  Felgar  (1963).  Specimens  of 
7075 -T6  aluminum  alloy  were  subjectsd  to  cyclic  loading  under  normal 
air  pressure  and  under  near-vacuum  conditions.  The  data  are  given  in 
Table  1. 

In  this  application,  near-vacuum  (simulating  outer  spa.  )  con¬ 
stitutes  the  use  environment  and  normal  air  pressure  the  ac  .  -lerated 
environment.  The  distribution  of  cycles-to-failure  in  vacu  :  is 
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Table  1 


ALLOY  FATIGUt  DATA 


Accelerated  Environment  (Air) 

Use  Environment  (Vacuum) 

Specimen 

10'  Cycles -to-Failure 

Specimen 

lO'1  Cycles-to-Failure 

X1 

1.23 

Y1 

3.3 

X2 

1.41 

Y2 

3.5 

X3 

1.G2 

Y3 

14.9 

X4 

2.08 

Y4 

40.7 

X5 

2.09 

Y5 

45.7 

X6 

2.30 

Y6 

46.8 

X7 

2.82 

Y7 

56.2 

denoted  by  G,  and  in  normal  atmosphere  by  F. 

In  our  example  the  number  of  observations  in  the  use  environment 

and  in  the  accelerated  environment  happen  to  be  the  same.  Often  there 

will  be  fewer  observations  available  in  the  use  environment. 

Using  formula  (3.2A  to  estimate  the  A  (X,)  's.  and,  mutatis  mutandis, 

n  1  - 

the  same  formula  to  estimate  the  v  tv  ^  >s  we  obtain 

m  i' 

X7(XX)  -  .1253 

*7(X2)  -  .2386 

A?(X3)  -  .3075 

A?(X4)  -  .4074 

A?(X5)  -  .5994 

A  AX.)  -  .o^mj 
7  o 

and 

Y?(Y1)  -  .0450,  1  -  I . b. 
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The  estimates  F,(t)  and  G,(t)  are  graphed  in  Figs.  1  and  2.  Since 
these  estimates  are  already  stochastically  ordered,  the  accelerated 
observations  give  us  no  further  information  concerning  the  distribution 
G.  A  statistical  test  for  IFRA ,  discussed  in  Sec.  4,  does  not  reject 
the  hypothesis  that  G  is  an  exponential  distribution.  This  is  also  made 

A 

credible  by  our  estimate  G^,  The  accelerated  data  on  the  other  hand  are 
significantly  IFRA  according  to  the  tota 1-time-on-test  statistic  dis¬ 
cussed  in  Sec.  4. 

3.4  ACCELERATED  LIFE  TESTS  ON  DUPLEX  CAPACITORS 

The  data  in  Table  2  were  obtained  from  conducting  accelerated  life 
tests  on  duplex  capacitors.  Acceleration  is  due  to  a  voltage  stress 
at  room  temperature.  The  criterion  for  failure  is  the  leakage  current's 
exceeding  a  certain  level.  The  times  to  failure  are  not  exact,  but 
approximately  close  to  the  actual  failure  times.  To  illustrate  our 
procedures,  we  consider  data  at  only  three  voltage  levels. 

Note  that  the  50-volt  test  was  censored  at  the  ninth  failure. 

Using  the  toUi-tirae-on-tes ;  statistic  described  in  Sec.  4,  the  50-volt 
and  70-volt  data  were  not  significantly  IFRA  relative  to  the  hypothesis 
of  exponentiality.  They  appear,  in  fact,  to  follow  an  exponential 
distribution.  The  60-volt  data,  on  the  other  hand,  are  significantly 
IFRA  even  at  the  .001  significance  level. 

As  we  would  expect,  the  data  Beem  to  be  stochastically  decreasing 
with  respect  to  Increasing  voltage.  A  statistical  test  for  stochastic 
ordering,  given  censored  data  and  using  a  modified  Wilcoxon  statistic, 


is  disc  «!8ed  in  the  next  section. 
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Table  2 

DUPLEX  CAPACITOR  ACCELERATED  LIFE  DATA 


50  Volts  (Room  Temp.) 

60  Volts  (Room  Temp.) 

70  Volts  (Room  Temp.) 

n  =  13 

n  =  13 

n  =  14 

Times  to  Failure  (hrs.) 

Times  to  Failure  (hrs.) 

Times  to  Failure  (hrs.) 

526 

56 

3 

1075 

517 

526 

1075 

517 

526 

1794 

517 

932 

2080 

517 

1023 

2891 

517 

1193 

2891 

517 

1223 

4211 

758 

2085 

4211 

758 

2296 

758 

2442 

758 

4216 

1016 

4216 

1186 

4216 

5589 

The  IFRA  estimates  of  the  survival  probability  function,  using 
Eq,  (3.4),  for  the  50- ,  60- ,  and  70-volt  data  are  graphed  in  Figs.  3, 
,  and  5,  respectively. 


4 


0  12  3 

t ,  thousands  of  hours 


Fig.  3  --  50-volt  data 


Survival  probability 
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STATISTICAL  TESTS  FOR  IFRA  (DFRA)  AND  STOCHASTIC  ORDERING 


Before  applying  the  techniques  of  this  paper  to  life  test  data, 
it  is  first  necessary  to  check  the  validity  of  the  assumptions  made. 

A  test  for  IFRA  (or  DFRA)  has  been  considered  by  Barlow  (1968)  and  is 
described  below.  (The  asymptotic  relative  efficiency  of  this  test 
against  various  alternatives  has  been  investigated  by  Bickel  and 
Doksum  (1967).)  Suppose  the  first  r  failures  out  of  a  sample  of  size 
n  are  recorded.  Let  S  ...  £  denote  the  observed  failure 

times.  Note  that  "withdrawals"  may  occur  between  X^  and  X^+^  and  that 
r  is,  in  general,  a  random  variable.  Let  n(u)  be  the  (random)  number 
of  items  on  test  at  time  u.  The  total  time  on  test  up  to  the  i  order 
statistic  is 

rxi 

T(X^)  *  j  n(u)du. 
r-1 

The  total-time-on-test  statistic  is  £  T(X,)/T(X  ).  Under  the  expo- 

i-1  1  r 


nential  hypothesis 


Z 


r-l  (  n 

l  T(X  )  -  T(X  ) 

1  1  * 

T(Xf)  /(r-l) / 12 


is  approximately  N(0,1)  even  for  relatively  small  r.  If  the  distribu¬ 
tion  of  time  to  failure  is  IFRA  (DFRA)  and  Y^  ^  Y£  <  •••  —  Yr  ^enotes 
an  independent  ordered  sample  from  an  exponential  distribution,  then 

r-l  r-l 

l  T(X  )  l  T(Y  ) 

_1 _ _  1  {<)  1-1 

T(X  )  st  st  T(Y  ) 


^Decreasing  failure  rate  average.  Cf.  Sec.  2. 
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where  a  denotes  stochastic  ordering.  Hence,  a  natural  test  rejects 
st 

exponentia lity  in  favor  of  IFRA  if 


r-1 

l  T(X  ) 

JL _ _ 

T(X  ) 


>  c 
—  a 


where  c  is  defined  by 
a 


r-1 

I  T(Y  ) 

JL _ _ 

T(Y  ) 


> 


c 

a 


a 


and  G(x)  -  1  -  e_x  for  x  >_  0. 

Table  3  presents  Z  values  corresponding  to  the  failure  data  in 
Tables  1  and  2. 


Table  3 

Z  VALUES  CALCULATED  FROM  THE  DATA  OF  TABLES  1  AND  2 


Data 

Z  Values 

Fatigue  Data  at 

air  pressure  (Table  1) 

2.8397 

Fatigue  Data  at 
vacuum  conditions 
(Table  1) 

0.4984 

Duplex  capacitors  at 

50  volts  (Table  2) 

0.7837 

Duplex  capacitors  at 

60  volts  (Table  2) 

3.2152 

Duplex  capacitors  at 

70  volts  (Table  2) 

0.5511 

The  conclusions  concerning  the  tests  of  exponentlclity  versus 
IFRA  for  the  alloy  data  (Sec.  3,3)  and  for  the  capacitor  data  (Sec.  3.4) 
were  reache  h-  comparing  the  Z  values  in  Table  3  with  percentage  points 
of  the  standard  normal  distribution. 
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TESTS  FOR  STOCHASTIC  ORDERING 

Although  stochastic  ordering  may  be  obvious  from  some  accelerated 
life  test  data,  there  may  be  instances  when  a  statistical  test  is 
appropriate.  For  example,  it  is  not  clear  that  there  is  the  desired 
stochastic  ordering  between  the  60-volt  duplex  capacitor  data  and  the 
70-volt  capacitor  data. 

A  generalized  Wilcoxon  test  for  comparing  arbitrarily  singly- 
censored  samples  has  been  proposed  by  E.  A.  Gehan  (1965  a,b)»  We 
assume  that  m  and  n  items  are  subject  to  life  test  and  we  observe 


X ' ,  ....  X*  r.  censored 

1  rx 

x  X  n-r.  failures 

r^+1  n  1 

Y.' ,  ....  Y1  r_  censored 

1  ’  r2  2 

Y  . ,  ,  ....  Y  m-r0  failures, 

tn  ^ 


(These  observations  are  not  necessarily  ordered.) 
Define 


/’I  <  Yj  or  X^  s  Yj 


1J 


0  if  X  -  Y  or  X'  <  Y  or  Y'  <  X  or  (X'  Y’)  ,  i.e.  ,  both 

*  ^  ^  are  censored 


+1  if  X,  >Y,  or  X'  2  Y 


and  calculate  the  statistic 


“  I  u 


where  the  sum  !e  over  all  run  comparisons.  To  test  the  hypothesis 


H 


o 


X  -  Y 
st 
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vei'Sus 


H,  :  X  s:  Y 
1  st 


we  reject  H  if 
J  o 


Z 


/fVar  (W  |  P.Hq)} 


c 

a 


where  is  determined  using  the  fact  that  Z  is  asymptotically  N(0,1). 


Var  (W 


P  *Hq) 


nm 


(n+m)  (n-hn-1) 


'i-1 


s 

l 

i=l 


‘iVV1’ 


1-1 


m^(m+n-M^-L^  p (n+m-3M^  ^ 


■1) 


where 


with  M  ®  0 
o 


with  L  *  0 
o 


and 

*  number  of  uncensored  observations  at  rank  i  in  rank  ordering 

of  uncensored  observations  with  distinct  values  (i  “  1,2,  s) ; 

*  number  of  right-censored  observations  with  values  greater  than 
observations  at  rank  i  but  less  than  observations  at  rank  (i+1). 

Var  (W  I  P ,H  )  is  the  variance  of  W  under  H  and  conditioned  on  the 
*0  o 


pattern  of  observations,  P. 
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5.  OPTIMALITY  OF  LEAST  SQUARE  IFRA  ESTIMATORS 

•J 

In  Sec,  3  we  described  a  method  for  estimating  an  IFRA  distribu¬ 
tion.  Since  Marshall  and  Prnschan  (n.d.)  have  shown  that  the  maximum 
likelihood  estimate  of  F  assuming  F  is  IFRA  is  not  consistent,  ve  have 
employed  a  least  squares  criterion  which  _i£  consistent.  Under  the  IFRA 
assumption  X(t)  *  is  nondecreas ing.  Given  the  observed 

failure  times  Xj  S  S  ...  S  based  on  a  random  sample  from  V  we  seek 
a  "good"  estimate  of  X(t).  Letting  F^  denote  the  product-limit  estimate 
(Eq.  (3. 1)) ,  then 


w 


-log  Fn(X1) 


provides  an  initial  estimate  of  X(X^)  which  will  enjoy  all  of  the 

properties  of  the  product-limit  estimate;  c.g.  ,  it  will  be  strongly 

consistent.  It  is  also  known  that  F  is  the  maximum  likelihood  estimate 

n 

of  F  in  the  class  of  all  distributions  (see  Kaplan  and  Meier  (lc-58)) 

and  hence  X  is  the  unrestricted  maximum  likelihood  estimate  for  X. 
n 

Define 


£  (X.)  =  Max  Min 


l  MVFJxji 


sSi  tSi 


l  FnfX^ 

j-s  n  j 


Clearly,  will  be  nondecreasing  at  order  statistics.  Furthermore,  if 
\  is  any  other  increasing  function,  then  \  has  the  property  thai 

„<X1>  -  ’-‘V'VV  j,  ‘VV  -  VV'VV  (5.1, 

i.=i  i=i 

+  [VX1>  *  >(X,)l2F„!Xi?- 
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Hence 


jj  ^‘V  '  S  t*‘V  - 


and  therefore  X^  is  closest  to  in  the  least  squares  sense  with  respect 
to  F^  in  the  class  of  nondecreasing  functions.  A  proof  of  inequality 
(5.1)  can  be  found  in  Brunk  (1965)  and  Marshall  and  Proschan  (1965, 
Theorem  5.1).  For  completeness  we  present  a  proof  for  the  discrete  case. 


Theorem  1:  (Marshall  and  Proschan) 

Let  h  (1  *  0,  1 . n)  be  a  nondecreasing  sequence,  an  ar¬ 

bitrary  sequence,  m  2  0,  i  *  0,  1 ,  2 ,  . . , ,  n  and 

t 

l  8  ,m, 

g  *  Max  Min  -  ; 

sSi  t2i  f 

l 

j-8  J 


?  (8,  ’  h  )2  m  >  l  (g.  -  h  ) 2  m  .  ♦  l  (g  -  g.)2  m^ 

l  1  i*0  i*0 


i.e.,  in  the  class  of  increasing  functions,  g  *  (g^,  ...,  gR)  is  closest 
to  g  *  (gQ,  gn)  in  the  least  squares  with  respect  to  the  measure 


Proof : 


Clearly,  g^  is  increasing  in  i.  We  need  only  show 


l  (4,  -  h  )<g  -  »  )  «  >  0. 
1*0 
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Suppose  is  constant  on  the  interval  [a,  b%  Then 

b 

I  (gj  -  hi)(g1  -  g£)  mi 
i-a 


Now 


b-1 

-  1 
i-a 

i 

l  (g,  -  S.)m 

j-a  3  J 

j  (hj+i  *  v  ■ 

b 

1  (gi  -  i>. 
J-a  3  3  J 

b-1 

\L , 

i  \ 

l  8jmj  / 

1-*.  .3-3.-  *  u 

i  8J  J 

v  r 

'L <hj+i  - 

hj)  + 

b 

I 

j 

b 

v 

( 

J-a  3  J 

i 

t 

l  »< 

a 

t 

l  8jmj 
j-a  J  J 
i 

l  «n1 

j-a  3 


l  8jmj 
>  Mini=#_ii 


t>i 


-1  l 

J-a 


m 


j 


Max  Min 
s^i  t>i 


l 

1^1 


J  J 


t 

l 

j-9 


m, 


*b 


*±- i  I 


Any  nondecreasing  function  defined  at  order  statistics  as  in  (3.2) 
is  obviously  permissible.  However,  for  definiteness,  we  define 


xn(t) 


0.  t  <  xx 

*n(xiJ *  xi  -  c  <  xi+l 


+  ®,  t  >  X  . 
—  n 


Note  that  thi9  estimate  will  be  closest  to  X  (t)  for  all  0  <  t  <  X 

n  —  n 

when  X  (X.)  -  X  (X.).  Our  estimate  for  the  distribution  function  will 
n  i  n  i 


then  be 
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i,  t  <  xx 

/> 

Fn(t)  -  expt-yx^t],  X1  fi  t  <  Xi+1 
0,  t  ?  X  . 

The  strong  consistency  of  F^  is  inherited  from  the  strong  consis¬ 
tency  of  F  .  By  Corollary  3.3  of  ”runk  (1965)  it  follows  that 
|X  (X  )  -  X (X±) }  £  e  for  i  -  1,2,  ....  n  implies  I ^n(X±)  -  A(X±) |  <  e 

for  i  =  1,2,  ...,  n.  Strong  consistency  follows  from  the  strong  con¬ 
sistency  of  \  . 

The  additional  inequalities  required  by  the  stochastic  ordering 
assumption,  G(x)  <_  F(x)  ,  are  immediately  evident  upon  examining  the 
graphs  of  -log  ^(t)  and  -log  F^Ct). 
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6.  IFRA  TIME  TRANSFORMATIONS 

j 

The  stochastic  ordering  assumption  on  the  time  transformation 
cr(t)  ^  t  is  probably  the  weakest  and  most  intuitive  requirement  on 
cr(t).  By  making  a  stronger  (but  more  difficult  to  justify)  assumption 
on  a(t)  we  should  obtain  "improved"  estimates. 

We  adopt  the  following  definition  of  an  IFRA  time  transformation:^ 

Definition 

<y(t)  is  an  IFRA  time  transformation  if  ar(t)/t  is  nondecre;sing  in 
t  5?  0  and  o(t)  5  t. 

This  definition  would  be  perhaps  quite  readily  acceptable  if  oft) 
depended  only  on  the  environment  and  not  on  the  original  dist:  button 
G.  Of  course  in  practice  this  seems  somewhat  unlikely.  Howevr,  this 
definition  leads  to  a  class  of  time  transformations  that  seems  in¬ 
tuitively  reasonable. 

Lemma  1. 

or(t)  is  an  IFRA  time  transformation  if  and  on  ly  if 
(i)  a(t)  is  nondecreasing  in  t  &  0; 

(ii)  a(t)  ?  t 

(ill)  for  every  IFRA  G,  we  have  that  F(t)  *  G[»(t)  is  IFRA. 

Proof: 

(a)  Suppose  o(t)  is  an  IFRA  time  transform,  tlon.  Then 
-log  F(t)/t  ■  [-log  G(o(t))/o'(t)][o'(t)/t]  is  none! ^creasing  in  !  i  0 
since  <y(t)  is  nondecreasing  in  t  2  0.  I.e.,  F  is  IFRA. 

^We  could  alternatively  consider  IFR  time  transformations.  A  dis¬ 
tribution  F  is  IFR  if  log[l  -  F(t)l  is  concave  on  its  support. 
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(b)  Suppose  conditions  (i)-(iii)  hold.  Choose  G  to  be  the  unit 

J 

exponential  distribution  and  define  F(t)  -  G[o(t)].  Then  o(t)/t  -• 
-log  F(t)/t  is  nondecreasing  in  t  s  0  since  F  is  I  FRA.  || 

Clearly,  if  o^(t)  and  02(6)  are  IFRA  time  transformations  then  so 
is  their  composition  o^Co^Ct)]. 

Lemma  2 

If  Y^  s  ^2  s  ***  s  ^n  ^Xl  s  ^2  S  **’  S  Xn^  are  or<*er  statistics 
from  G  (F)  and  F(t)  *  G[or(t)],  where  o  is  an  IFRA  time  transformation, 
then 

st 

Yj  "  Yi  -  Xj  -  X,  *>r  i  <  J 

at 

where  >_  means  stochastically  greater  than.  In  particular,  the  sample 
range  from  G  is  stochastically  greater  than  the  sample  range  from  F. 


Proof: 

s  t 

Let  Y^  ■  G  ^F(X^)  and  note  (Y^,  ...,  Y^)  ■  (Y^, 

_  1 

x  c/  *■  ^  a(t) 

- 1”1”  "  — C  —  ^  ant^  nondec teasing  in  t  >_  0, 


G_1F(Xj)  -  G"1F(Xi) 
Xj  “  Xi 


Y  ) .  Since 
n 


which  implies  the  result.  | | 
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From  the  results  of  Barlow  and  Pro3chan  (1966),  it  is  easy  to 

•j 

verify  that  IFRA  time  transformations  stochastically  decrease  order 

statistics  and  the  sample  coefficient  of  variation.  It  follows  that 

population  coefficient  of  variation  and  moments  are  also  decreased. 

Assuming  that  a(t)/t  **  G  ^F(t)/t  is  nondecreasing  in  t  ^  0  is 

equivalent  to  assuming  that  F(t)  crosses  G(0t)  at  most  once  and  from 

below  if  at  all  for  every  0  >  0.  In  this  sense,  F  is  "sharper"  than  G. 

In  analogy  with  the  least  squares  estimates,  we  confine  attention 

to  estimates  of  the  form  G  and  F  where 

v  Y 


-log  G^ (t) 


0  ,  0  <  t  <  Y1 

£Y(Y1),  <  t  <  Y1+1 


and  -log  F^(t)  is  defined  similarly. 

Given  observations  Y  <  Y_  <  . . .  <  Y  (X,  <  X„<  . . .  <  X  )  from 

1—2—  —  m  1—2—  —  n 

a  distribution  G  (F)  we  wish  to  estimate  both  F  and  G  assuming 

(i)  F  and  G  are  IFRA, 

(II)  G  *F(t)/t  is  nondecreasing  in  t  >_  0 
(iii)  G-1F(t)  >_  t  for  t  ^  0  . 

* 

Define  y(t)  -  -log  G(t)/t  and  y  (t)  -  -log  G  (t)/t  where  G  is 

mm  m 

the  usual  empirical  distribution  or  product  limit  based  upon  observa¬ 
tions  Y,  <  Y-  <  ...  <  Y  from  the  normal  or  "use"  environment.  We 
consider  the  problem  of  minimizing 
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o' 

l 


lY(x)  -  Y_(x) ]  dG  (x) 
m  m 


Or> 

+/„ 


[X(x)  -  X(x) ]  dF  (x) 
n  n 


(6.1) 


with  respect  to  y  and  X  subject  to 


0  <  y(Vx)  1Y(Y2)  1  ...  <  Y(Yn), 

0  1  X(XX)  <_  X(X2)  1  ...  1  X(Xn)  , 

G  *F.(t)/t  is  nondecreasing  in  t  2  0,  and  G”*F.(t)  S  t.  If  G~*F  (t)  is 
y  \  Y  \  n  n 

nondecreasing  in  t  >  0  and  G”*Fn(t)  ^  t  then  clearly  the  least  squares 

solutions  y  and  X  also  minimize  (6.1). 
tn  n 

To  understand  the  restrictions  on  the  time  transformation,  we 
make  use  of  the  so-called  Q  -  Q  plot  (Q  for  quantile)  which  is  widely 
used  by  data  analysts  [cf.  Wilk  and  Gnanadesikan  (1968)].  Let  G^  (F^) 
denote  estimates  of  G  (F) .  A  Q  -  Q  plot  based  on  G  and  F1  is  merely 
a  plot  of  G~^F^.  If  G^  and  F^  are  the  empirical  distribution  functions, 
this  provides  a  quick,  although  heuristic,  check  on  our  assumptions 

_  i  .v  —  2. 

concerning  G  F.  Let  H(x)  »  1-e  and  consider  the  graph  of  H  G^ 

and  H  ^F^  as  illustrated  in  Figs.  6  and  7.  Note  that 

[H_1G  ]"lH"1F.  -  G“1F. . 
y  A  y  A 

The  parameters  associated  with  line  segments  are  their  slopes. 

The  Line  segments,  if  extended,  would  pass  through  the  origin.  Fix  an 
ordinate  value,  say  v.  Define  x(v)  to  be  the  largest  value  of  x  such 

that  H~*F^(x)  ^  v  and  y(v)  to  be  the  largest  value  of  y  such  that 

H  '(My)  v.  Letting  v  range  over  all  positive  values,  plot  'he 
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locus  of  points  (x(v) ,y(v))  in  Fig.  8.  Since  G_1F(t)/t  is  nondecreasing 

in  t  >  0,  we  demand  the  same  of  G  ^F.  (t)/t.  Hence  G  F.  will  consist 

—  Y  A  Y  A 

of  line  segments  with  increasing  slopes,  each  greater  than  or  equal  to 

one.  Since  G_1F,  cannot  be  flat  over  x-axis  intervals,  we  must  have 

Y  * 


X1X!  i*lY! 


(6.2) 


If  we  were  to  plot  points  y^Y^,  (i  *  1,2,  . . . ,  m)  and 

X  jXj,  ^jXj»  <J  *  1,  2,  ...»  n)  on  the  same  line,  we  would  see  that  in¬ 
tervals  of  the  form  (y^  jY^ ,  YjY^)  must  fall  within  intervals  of  the  form 
(X  xXj,  XjXj)  in  order  that  g"^  not  be  constant  over  intervals. 

To  minimize  (6.1)  subject  to  G,  F  I  FRA  and  G_1F(t)/t  nondecreasing 
in  t  >  0  we  must  minimize  (6.1)  subject  to  the  constraints 


where  we  must  consider  every  possible  positioning  sequence  (j^tjj . Jm) 

in  turn.  Given  a  positioning  sequence  ( J  ^ » J  2 »  •••»  Jm)  we  must  aJso 
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A  computer  program  Is  necessary  to  solve  this  problem  in  general.  If 
n  Is  large  so  that  we  have  a  good  estimate  of  F,  then  we  may  let 
\j  “  Aj  and  .minimise,  (6.1)  with  respect  to  y  subject  to  the  restrictions 
on  y . 

It  is  intuitively  clear  from  Fig.  8  that  if  the  y^  estimates  are 
too  high  the  slopes  A^/y^  may  not  satisfy  our  requirement  A^/y^  ■  1. 
Hence  the  A^'s  will  tend  to  force  our  estimates  of  the  y's  down. 
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7.  A  MODIFIED  PROBLEM 


By  introducing  additional  restrictions  on  G,  the  failure  distribu¬ 
tion  for  the  use  environment,  we  may  obtain  an  explicit  solution  to 
•■he  problem  of  Sec,  6,  Assume 


G(t) 


t  <  Y, 


b  exp[ - v( t  -  Y^  ,  t  >  Yj 


where  b  and  v  are  parameters  to  be  determined  subject  to  the  restriction 
e ~^1  s  b  £  1.  The  restriction  on  b  insures  that  G  is  IFRA.  This  is  a 
more  stringent  assumption  on  G  than  made  previously.  However,  one  can 
imagine  situations  where  it  might  be  approximately  true.  Figure  7  is  now 
a  single  ray  starting  above  Y^  with  slope  V  and  continuing  indefinitely 
to  the  right. 

We  consider  the  problem  of  minimizing 


I 


(Y  - 


Y_(x) ]' 
m 


dC  (x) 
m 


or 

X 


[X(x)  - 


Xn(x)]‘ 


dFn(x) 


(7.1) 


with  respect  to  y  and  X  subject  to 


0  <  y 


Xlh 


~  Y, 


Note  that  in  most  practical 
will  automatically  imply  1 


a  ^  n 

and  1  <  —  <  — 
-  Y  “  Y 

situations  Y 
X1 

<  —  .  Let  6 
Y  o 


-  <  ...  < - . 

-  -  Y 

>1*1 

>  X  so  that  y  s  — 
Y1 

=  Yxy  and  6j  -  X^, 


j-1,2,  . . .  ,  n  -  1 . 

We  now  assume  that  we  arc  dealing  with  a  complete,  noncensored 
sample  with  no  tied  observations  and  no  losses.  Rewriting  (7.1)  we 


33- 


subject  to  0  <  6,  <  ...  <9  , .  The  solution  is  immediate  from 

J  o  —  1  —  —  n-1 

Theorem  1  if  we  make  the  following  identification: 

\  '»(V 

8o  "  Y1  Im  5  8j  *  XlVX]}  for  J  *  2»  "-1 

and 

1  1  c  .  ,  „ 
m  ■  T~,  ■ — 7  for  j  »  1,  2,  n. 

n  J  nXj 

Applying  the  method  of  this  Section  to  the  fatigue  data  of  Table 
1  we  find 

\  -  .0193 

\l  «  .1253 

X2  *  .2386 

-  .3075 

X.  «  .407 4 
4 

X5  -  .5994 
\6  *  .8460 
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Comparing  these  estimates  with  those  in  Sec.  3.3,  we  note  that  our 

■j 

estimate  of  the  failure  rate  in  vacuum  (the  use  environment  in  this 
case)  is  less,  i.e.,  y  ■  .0193  versus  .0450  in  Sec.  3.3. 
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